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ABSTRACT 

We measure the angular clustering of ^ 2000 radio sources in The Bootes Deep Field, cover- 
ing 5.3 deg^ down to S'i.4ghz = 0.2 mJy. With reference to work by Blake & Wall, we show 
that the size distribution of multi-component radio galaxies dominates the overall clustering 
signal, and that its amplitude extrapolates smoothly from their measurements above 5 mJy. 
The upper limits on any true galaxy-galaxy clustering are consistent with the clustering of 
sub- mJy radio-loud AGN being effectively diluted by the more weakly-clustered IRAS-type 
starburst galaxies. Source count models imply that the survey contains ~ 400 of the latter 
galaxies above 0.2 mJy out to z ^ 1 — 2. Measurement of their clustering must await their 
identification via the optical and infrared data due on this field. 
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1 INTRODUCTION 

Measurements of the clustering of galaxies can shed light on the 
formation of large scale structure and of the galaxies within it. In 
the local Universe (z < 0.1), optically-selected galaxies are un- 
biased tracers of the mass, following a spatial correlation func- 
tion of the form ^(r) — (r/ro)^^, with ro — 5Ah^^ Mpc and 
7 = 1.8 (Davis & Peebles 1983). Results from the 2dF QSO red- 
shift survey show that al z — 1.3, QSOs are comparably clus- 
tered (ro = 5.7h~^ Mpc; Shanks et al. 2001). In contrast, the 
local population of IRAS starburst galaxies are less strongly clus- 
tered (ro — 3.8h~^ Mpc; Saunders et al. 1992), whilst the op- 
posite is true of early-type galaxies (for local L > L* ellipticals 
ro = 9 - llh'^ Mpc, e.g. Guzzo et al. 1997, Willmer et al. 1998). 

Extending such work to higher redshift populations can con- 
strain the galaxy systems into which they evolve: e.g. Adelberger 
(2000) measured the clustering of ^ ~ 1 Balmer-break galaxies and 
derived ro = S.Oh^^ Mpc (comoving). This is consistent with the 

2 ~ 1 Balmer-break galaxies being both unbiased tracers of mass 
and the progenitors of local L* galaxies. The comparable comoving 
space densities of the two populations supports this interpretation. 
In contrast, at z ~ 3 the Lyman-break galaxies are more strongly 
clustered (see Porciani & Giavalisco 2002), implying that they are 
strongly biased relative to the mass distribution. Similarly, analy- 
sis of the extremely red objects (EROs) demonstrates that passively 
evolving early-type galaxies at z ~ 1.2 have ro — 12±3/i^^ Mpc, 
comparable to that of local L* ellipticals (Daddi et al. 2001; Mc- 
Carthy et al. 2001). This rules out a scenario in which Lyman-break 
galaxies evolve into local bright ellipticals, and requires that the 
bias evolves with redshift such that the comoving density of el- 
lipticals is invariant to z ^ 1 and beyond (see also Moustakas & 
Somerville 2002). 

In contrast, radio surveys probe galaxy clustering over larger 



scales, as they cover large areas of sky with a broad redshift dis- 
tribution. This does, however, reduce the clustering signal when 
projected on the sky. Nevertheless, numerous positive measure- 
ments of the angular clustering of radio sources have been made. 
The most complete analysis to date is that of Blake & Wall (2002) 
(hereafter BW02), who measured the angular correlation function 
of sources in the 1.4 GHz NVSS from flux thresholds of 50 down to 
5 mJy (see also Overzier et al. 2002). They found that the clustering 
at small angles is due to the size distribution of multi-component 
radio galaxies, and that the galaxy-galaxy correlation function has 
the universal slope of 7 = 1.8, and (under certain assumptions) a 
correlation length of ro ~ 6h^^ Mpc. 

Below ~ 1 mJy (at 1.4 GHz) there is evidence for the appear- 
ance of a different radio source population which may replace the 
radio-loud AGN as the dominant population below a few 100/iJy 
(Windhorst 1984, Windhorst et al. 1985, Benn et al. 1993). These 
sources are more distant analogues of the dusty starburst galaxies 
selected by IRAS, and models of the radio source counts require 
that this population undergo substantial luminosity evolution out 
to z ~ 1 (Rowan-Robinson et al. 1993; hereafter RR93). Mea- 
surements of the clustering of the sub-mJy population could thus 
extend the local IRAS clustering measurements of starburst galax- 
ies out to jz ~ 1, for comparison with those of the Balmer-break 
galaxies, the EROs referred to above and the dusty, star-forming 
EROs at 2 ~ 1 which have ro no larger than 2.5h^^ Mpc (Daddi 
et al. 2002). Such a measurement was attempted by Georgakakis et 
al. (2000) for sources with Si. 4 > 0.5 mJy in the ~ 3deg^ Phoenix 
radio survey (Hopkins et al. 1998). Whilst they found some evi- 
dence for clustering, the uncertainties were such that they could 
not determine whether the starbursts have a lower ro than the radio- 
loud AGN. 

To overcome these limitations we have initiated The Bootes 
Deep Field, a radio survey covering ~6 deg^ down to a limiting 
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5(7 sensitivity of 140 /iJy at 1.4 GHz (de Vries et al. 2002). The 
source finding algorithm identifies 3172 distinct sources, of which 
~ 10 per cent are resolved by the 13 x 27 arcsec beam (the actual 
number of sources used in our angular clustering analysis is smaller 
than this, due to our use of a subset of the data with a slightly 
higher flux limit, and the use of certain algorithms to identify mul- 
tiple component sources; see section 2). It will be complemented 
by imaging in six optical and near-infrared wavebands (as part of 
the NOAO Deep Wide-Field Sur-vey), and at longer infrared wave- 
lengths with SIRTF. These data will enable morphological classifi- 
cation of the radio sources into starburst or AGN, photometric red- 
shift estimation, and an examination of the far-infrared:radio corre- 
lation for starburst galaxies out to cosmological distances (building 
on work by Garrett 2002). Direct measurement of the real space 
correlation functions of both radio-loud AGN and starbursts (and of 
the cross-correlation between them) will also be possible to 2; ~ 1. 
Furthermore, these measurements will be free from the confusion 
caused by multi-component radio sources in existing catalogues 
without identifications. Until this dataset is assembled, however, 
we are confined to use of the angular correlation function, and it is 
these measurements and their interpretation that we present here. 



2 ANGULAR CLUSTERING MEASUREMENTS 

Full details of the Bootes survey characteristics and source ex- 
traction procedure are given in de Vries et al. (2002). As detailed 
therein, the survey is complete down to a limiting flux density of 
0.2 mjy, as determined by the point at which the source counts 
exhibit the first systematic deviation from a low-order polynomial 
fit (see their Fig. 9; N.B. The source counts shown in de Vries et 
al. are all too low by a factor 2.035 due to a numerical error in the 
construction of this plot). Since extracted sources must have peak 
fluxes at least 5 times higher than the local rms noise level, for the 
clustering analysis we exclude the comers of the survey area where 
the latter exceeds 0.04 mJy; the remaining area covers 5.32 deg^. 

The angular two-point correlation function is defined 

through the expression for the probability, 5P, of finding two 
sources in solid angle elements Sfli and 5Q2, with angular sep- 
aration 6: 



SP = N'^SQiSQiil + uj{0)), 



(1) 



where A*' is the mean areal source density. The function uj{0) 
is therefore a measure of the deviation from a random distribution. 
It can be calculated by comparing the number of source pairs within 
a given range of angular separation with the number of pairs in a 
large random catalogue covering the same area. Of the variety of 
estimators for u!{0) which have been proposed, we follow Ren- 
gelink (1999) and use that due to Hamilton (1993) because of its 
robustness at large angular scales: 



DD ■ RR 



(nu - l)(nR - 1) DR-DR 



(2) 



where DD is the number of data-data pairs within the angu- 
lar bin centred on separation 9, and RR and DR are the numbers 
of random-random and data-random pairs, respectively, within the 
same separation interval. We use a random catalogue containing 
riR = 25000 sources, vastly exceeding the size of the data cat- 
alogue (no)- Errors on the individual w{8) points are computed 
with the method of Ling, Barrow & Frenk (1986), by calculating 
the standard deviation in uj{6) among 20 pseudo-random resamples 



of the observational dataset. Such estimates exceed the Poisson er- 
rors, which are correct only for unclustered data. 

Using this procedure, u}{d) was computed in equally-spaced 
logarithmic angular intervals between 0.25 and 155 arcmin for sev- 
eral subsamples with lower limiting fluxes between 0.2 and 2 mJy. 
Each was then fitted with a function of the form: 



a;(6»)flt = Ae-" - C 



(3) 



with 6 in degrees. The quantity C (the integral constraint; 
Groth & Peebles 1977) is a bias resulting from the finite bound- 
ary of the survey and is given by: 



uj{e)dQ,idQ.2. 



(4) 



For the conventional power-law form ijj{0) — A6^^, Monte- 
Carlo computation of the integral yields C — 1.1 54 A and 1.392 A 
for 5 = 0.8 and 1.1, respectively. 

A complication in the analysis of radio source clustering is 
the signal at small angular separations (below a few arcminutes) 
caused by sources with multiple components. This problem will 
be largely alleviated when such sources can be reliably identified 
with the optical and infrared data. In the meantime, we begin by 
adopting some prescriptions used by Georgakakis et al. (2000), and 
before that by Magliocchetti et al. (1998) and Cress et al. (1996), 
to identify the genuine multi-component sources prior to measur- 
ing the angular correlation function. Using the 6 oc relation 
found by Oort (1987), we consider as a single object all doubles 
with 6 < 20\/Ftotai, where 6 is their separation in arcseconds 
and Ftotai their summed flux density in mJy; furthermore, we only 
collapse doubles whose individual component fluxes differ by less 
than a factor of 4 (since components of genuine doubles are ex- 
pected to have correlated fluxes; Magliocchetti et al. 1998). Of the 
214 double sources identified by the source extraction algorithm, 
50 are identified as genuine by these criteria. The 48 sources with 
3 or more components (39 triples, 9 quadruples) were examined 
by eye and subjective criteria (e.g. relative component fluxes, mor- 
phology) were used to decide whether or not to treat them as a 
single source, resulting in the assignment of 78 separate sources. 

Fits to the data were performed between separations of 1.5 and 
20 arcmins, and are shown in Fig. [l] for the 0.2 and 2 mJy subsam- 
ples. The fitted parameters for these and other subsamples are listed 
in Table 1 . As expected, the amplitudes are lower when fitted with 
steeper power-laws. The variation of amplitude with flux limit for 
our survey is shown in Fig. ^ Although the error bars are large, our 
measured amplitudes of ~ 0.01 for flux limits of 1-2 mJy are con- 
sistent with the results of Georgakakis et al (2000) for the Phoenix 
survey, and with those of Cress et al. (1996) who measured an am- 
plitude of ~ 0.008 for FIRST sources in the flux interval 1-2 mJy. 

However, clustering analyses of the NVSS and FIRST radio 
surveys by BW02 and Overzier et al. (2002) are at odds with the 
earlier measurements from the FIRST survey. Both papers ques- 
tion the efficacy of the procedures previously used (which we also 
adopt) to identify the multi-component sources: they find that the 
angular correlation function below 6 arcmin is dominated by multi- 
component sources and that the true cosmological clustering am- 
plitude is essentially constant at ~ 10"^ from 3 mJy to ~ 50 mJy. 
It seems implausible that the amplitude could drop by a factor of 
10 from flux limits of 1-2 mJy to 3 mJy. We thus conclude that the 
increase in amplitude from 0.2 to 2 mJy in Fig. ^ is due to the resid- 
ual effects of multi-component sources at the higher fluxes. Indeed, 
Fig. |l| shows excess signal above the fitted power-law at small an- 
gles for the 2 mJy subsample. 
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Table 1. Fitted amplitudes of the angular con'elation function after at- 
tempted removal of multi-component sources. The errors are Icr. 



Flux density ( mjy) 


10^yl(<5 = 0.8) 


10^A(5 = 1.1) 


> 0.2 


2.1 ± 1.6 


0.86 ± 0.67 


> 0.4 


2.0 ±3.0 


0.86 ± 1.3 


> 0.6 


4.5 ± 3.9 


1.9 ± 1.6 


> 0.8 


4.2 ± 5.4 


1.8 ±2.3 


> 1.0 


4.9 ± 4.8 


2.3 ± 2.0 


> 1.5 


9.4 ± 7.7 


3.9 ±3.2 


> 2.0 


12 ± 8.3 


5.4 ±3.6 



Flux density > 0.2 mJy 




Flux density > 2.0 mJy 




10 

(9(arcmin) 



— 



100 



Figure 1. Angular correlation functions for subsamples with lower 
flux limits of 0.2 and 2 mJy, after the attempted removal of genuine 
multi-component radio sources, as described in section 2. The sohd 
lines show fits to the data over 1.5 < S < 20 arcmin with the 
function given in eqn. 3, for 5 = 0.8. There is an excess of data 
over the model at the smallest separations in the 2 mJy data, likely 
due to the residual effects of the multi-component sources. 
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Figure 2. Variation of the fitted amplitude of the angular correlation 
function with flux limit in The Bootes Deep Field (after attempted 
removal of multi-component sources from the catalogue), for as- 
sumed power-law indices of 5 = 0.8 (filled circles) and 5 — 1.1 
(open squares); for clarity the latter points have been displaced 
slightly in flux. 



Table 2. Fitted amplitudes of the size-distribution and cosmological clus- 
tering power-laws derived from a two component fit. The errors are Icr. 



Flux density 
(mJy) 


size power- 
10'^ A{5 = 


law 
3.4) 


cosmological clustering 
lO'^i A(S = 0.8) 


> 0.2 


2.4 ± 1.2 




0.55 ± 1.4 


> 0.3 


1.9 ± 1.6 




0.38 ± 1.9 


> 0.4 


2.2 ± 1.7 




1.2 ±2.5 


> 0.6 


4.1 ± 2.1 




1.5 ±3.7 


> 0.8 


3.4 ±3.4 




2.1 ±4.4 


> 1.0 


4.8 ± 5.3 




2.5 ± 4.7 


> 1.5 


9.7 ± 7.6 




4.8 ±6.9 


> 2.0 


14 ±8.3 




4.3 ± 8.5 



In an attempt to overcome the effects of multi-component 
sources, we instead apply none of the above component combining 
procedures and fit the resulting correlation functions with the sum 
of two power-laws of the form shown in eqn (3). The first has 5 — 
3.4 (fixed) and is due to the multi-component sources which dom- 
inate the signal below 6 arcmin, as found by BW02 and Overzier 
et al. (the corresponding integral constraint is A = 1191C, with 
the correlation function truncated below the angular resolution of 
the survey to keep the integral in eqn. 4 finite). The second has 
5 = 0.8 (fixed) and is due to the cosmological clustering. Fits are 
performed to data above 1 arcmin, as shown in Fig. ^ for the 0.2 and 
2 mJy sub- samples. Table 2 lists the results. There is now no sig- 
nificant detection of cosmological clustering and the effects of the 
multi-component sources dominate the w{8) signal. Furthermore, 
Fig. ^ shows that the amplitude of the size distribution power-law 
is consistent with the extrapolation to fainter fluxes of the 1/a de- 
pendence (a being the surface density of radio sources) found from 
5 to 50 mJy by BW02. With reference to their section 3, this is con- 
sistent with the quantity e/n varying by no more than a factor of 
~ 2 from 50 to 0.2 mJy, as judged from the scatter around the fit 
in Fig. ^ (e being the fraction of sources observed to have mul- 
tiple components, n the average number of radio components per 
source). 
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Figure 3. Angular correlation functions for subsamples with lower 
flux limits of 0.2 and 2 mJy, without applying any component com- 
bining procedures to the catalogue; the solid lines are fits to the 
data above 1 arcmin with a double power-law, as described in the 
the text. Note the different y-axis scales. 



3 INTERPRETATION 

We now model the evolution of the space density and clustering 
of the sub-mjy radio source population, to determine in par- 
ticular the effect on the clustering signal of the appearance of 
the starburst galaxies. The spatial and angular correlation func- 
tions, ^(r) and u-i{6), are related via the relativistic Limber equa- 
tion (see e.g. Peebles 1980; Magliocchetti et al. 1999); when 
^(r, z) = (r/ro)~''(l + z)~^^'^^^ (r and ro being proper lengths), 
w{6) = Ad~^''~^'^ {9 in radians) where for a spatially flat cosmol- 
ogy with non-zero cosmological constant (SIa + = 1): 



7-3- 



''Q{z)dz 



with 



B 



{irN{z)dz) 



r(V) ' 



Q{z) = [{l + zf + n^^-l] 



(5) 

(6) 
(7) 




2000 



Figure 4. The circles represent the amplitude of the contribution to 
Lo{8) from the multi-component sources, from Table 2. is the 
number of sources above the corresponding flux limit (oc a, the 
surface density - with no component-combining procedures ap- 
plied). The fit to these points with Adoubie oc A'^" is shown (with 
a = —0.99 ± 0.6), and is consistent with the A'^"^ extrapola- 
tion from the higher flux limits found by BW02 (whose point for 
sources > 5 mJy is indicated by the star). 



/o 



(8) 



and N(yz) is the redshift distribution of sources above the flux 
limit. When a sample comprises two sub-populations, A and B (in 
this case AGN and starbursts), with different clustering properties, 
the signal for the whole is given by: 



- 2fAfBl^AB, 



(9) 



where /a and /b are the fractions of the population in the two 
classes, and uja and lub their individual correlation functions. The 
term ujab is the cross-correlation between the two populations. In 
our analysis, we use the Limber equation to compute separate io{9) 
for the AGN and starbursts, and then combine them using eqn. (9). 
For simplicity, we set the cross-correlation to zero. 

Redshift distributions for each flux-limited subsample were 
computed using the Dunlop and Peacock (1990) pure luminos- 
ity evolution model for the combined population of flat (q = 0; 
Si, oc v") and steep (a = —0.8) spectrum AGN (taking the pa- 
rameters from their table C3, shifted to 1.4 GHz). For the star- 
forming galaxies, we use the determination of their local 1.4 GHz 
luminosity function given by Sadler et al. (2002), with pure lumi- 
nosity evolution of the form (1 + z^^ out to 2 = Zcut, with no 
further evolution thereafter (RR93); with Q = 3.1 (as in RR93), 
we find that ^cut = 1.5 can reproduce the integral source counts 
satisfactorily, as shown in Fig. n (we plot integral, rather than dif- 
ferential, source counts as the former more clearly show the num- 
ber of sources within each flux-limited sample used for the clus- 
tering analysis). For flux thresholds varying from 0.2 to 2.0 mJy, 
the model predicts that the starburst fraction in the integral source 
counts decreases from around 0.2 to 0.03. Fig. |^ shows the model 
redshift distributions for a flux limit of 0.2 mJy. It should be noted 
that the source counts (in Fig. ^ and redshift distribution (in Fig. |^ 
and used as an input to eqn. 5), were obtained with a cosmology 
of Ho — 50kms^^ Mpc~^ and Q,m = 1.0, for which the above 
luminosity functions were also derived. However, redshift distribu- 
tions and source counts are observational quantities, and thus inde- 
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0.2 0.5 

Flux limit (mJy) 



Figure 5. The points show the integral source counts in the central 
5.32 deg^ of The Bootes Deep Field where the catalogue is com- 
plete down to 0.2 mJy. The dashed and dotted lines show the con- 
tributions of the AGN and the starbursts, respectively, and the solid 
line their sum. The ^/N error bars are smaller than the symbols. 
See the text for details. 



Sj ^ > 0.2 mJy 




Figure 6. Model redshift distribution for flux limit of 0.2 mJy. The 
solid and dashed lines denote the starbursts and AGN, respectively. 



pendent of the assumed cosmology. This is not, therefore, inconsis- 
tent with our use of a different cosmology (S7m = 0.3, JIa = 0.7) 
in the clustering calculations which follow. 

Using the above formalism, we compute the amplitude of the 
total angular correlation function for a baseline model in which 
ro = 6.0/i~^ Mpc for the AGN (comparable to that measured by 
BW02), and ro = 3.0/i^^ Mpc for the starbursts (within the range 
measured for local IRAS starburst galaxies). For both correlation 
functions, we take 7 = 1.8 and the clustering evolution param- 
eter e = 7 — 8; the latter corresponds to constant clustering in 
co-moving coordinates, as observed for elliptical galaxies (see sec- 
tion 1). We take Q,m = 0.3 and f^A = 0.7 and note that the re- 
sults are independent of Hq. In Fig. ^ we show, as functions of the 
flux limit, the amplitudes of the angular correlation functions of 
the AGN and starbursts separately, and for the combined popula- 
tion (combining the two signals as in eqn. 9). The upper limits in 
this figure are our 90 per cent confidence limits on the amplitude 
of the galaxy-galaxy clustering, obtained from a two power-law fit 
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Figure 7. The solid line shows the model-predicted amplitude of 
the angular correlation function as a function of flux limit. The 
predicted contributions of the AGN and starbursts alone (i.e. not 
weighted by the squares of the population fractions as shown in 
eqn. 9), are denoted by the dashed and dot-dashed lines, respec- 
tively. The arrows denote the 90 per cent confidence upper limits 
on the galaxy-galaxy clustering amplitude; the star is the measure- 
ment from Overzier et al. (2002) using FIRST and NVSS data. 



with the amplitude of the size distribution power-law fixed at the 
value set by the 1/a extrapolation of the 5 mJy point of B W02 ( see 
Fig. The constancy of the model AGN clustering amplitude re- 
flects the invariant shape of the AGN redshift distribution over this 
range in flux limit. The large amplitude of the starburst clustering 
at the high flux end reflects the fact that the brightest starbursts 
are local objects with a fairly narrow redshift distribution; moving 
towards fainter fluxes the redshift distribution broadens, reducing 
the amplitude. Given the small starburst fraction in the model over 
this flux range, the effective signal for the combined population is, 
with reference to eqn. (9), essentially ujcs — ./agn'^agn- There is 
thus little scope for measuring the clustering of starburst galaxies 
alone with our present data. In order to do so using unidentified ra- 
dio data one would have to push down well below O.I mJy to the 
level at which the starburst fraction becomes substantially higher. 
However, when the ~ 400 IRAS-type starburst galaxies above the 
0.2 mJy flux limit (see Fig. ^) can be identified in the follow-up 
data due on this field, measurement of their clustering strength will 
become easier. 



4 CONCLUSIONS 

The Bootes Deep Field is a 2.5 x 2.5 deg region which has been 
surveyed at 1.4 GHz down to an rms noise level of 28 /ijy at its 
centre. The area will be covered at 325 MHz, in six optical and 
near-infrared bands, and at longer infrared wavelengths as part of 
a SIRTF legacy programme. Here we have measured the angular 
two point correlation function of the 1 .4 GHz sources down to the 
survey limit of 0.2 mJy. We find that the size distribution of multi- 
component radio galaxies dominates the overall signal even at these 
faint fluxes, with an amplitude consistent with the extrapolation of 
the 1/(7 variation established over flux limits from 5 to 50 mJy by 
BW02 (cr being the surface density of radio sources). This implies 
that the fraction of multi-component sources (normalised by the 
average number of components per source) varies by no more than 
a factor of 2 from 50 mJy to 0.2 mJy. 
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Only upper limits can be placed on the strength of any true 
galaxy-galaxy clustering. These limits are, however, consistent 
with the extrapolation of the amplitudes measured at higher fluxes 
by e.g. BW02, and with a model in which the clustering of radio- 
loud AGN is effectively 'diluted' by the more weakly clustered 
starburst galaxies. Source count models imply the latter popula- 
tion comprises 20 per cent of the sources above 0.2 mJy (some 400 
objects), with a broad redshift distribution peaking at z ~ 1 — 2. 
Measurement of the clustering of these galaxies alone would ex- 
tend earlier IRAS measurements from z ^ 0.2, but must wait until 
they can be identified with the follow-up data on this field. 
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